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Abstract 

We employ Landau's theory of normal Fermi liquids to study the quasiparticle interaction in nuclear 
matter in the vicinity of saturation density. Realistic low-momentum nucleon-nucleon interactions 
evolved from the Idaho N 3 LO chiral two-body potential are used as input potentials. We derive for 

(N 

the first time exact results for the central part of the quasiparticle interaction computed to second 

3 . 

order in perturbation theory, from which we extract the L = and L = 1 Landau parameters as 

00 

well as some relevant bulk equilibrium properties of nuclear matter. The accuracy of the intricate 
numerical calculations is tested with analytical results derived for scalar-isoscalar boson exchange and 
(modified) pion exchange at second order. The explicit dependence of the Fermi liquid parameters on 
the low-momentum cutoff scale is studied, which provides important insight into the scale variation of 

C ' 

phase-shift equivalent two-body potentials. This leads naturally to explore the role that three-nucleon 
forces must play in the effective interaction between two quasiparticles. 
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I. INTRODUCTION 



Describing the properties of infinite nuclear matter has long been an important benchmark 
for realistic models of the nuclear force and the applied many-body methods. Recent calcula- 
tions jl-4] have shown that the (Goldstone) linked-diagram expansion (up to at least second 
order) can provide an adequate description of the zero-temperature equation of state when real- 
istic two-nucleon and three-nucleon forces are employed. In the present work we study nuclear 
matter from the perspective of Landau's Fermi liquid theory js - 8 1 , which is a framework for de- 
scribing excitations of strongly-interacting normal Fermi systems in terms of weakly-interacting 
quasiparticles. Although the complete description of the interacting many-body ground state 
lies beyond the scope of this theory, various bulk equilibrium and transport properties are 
accessible through the quasiparticle interaction. 

The interaction between two quasiparticles can be obtained microscopically within many- 
body perturbation theory by functionally differentiating the total energy density twice with 
respect to the quasiparticle distribution function. Most previous studies using realistic nu- 
clear forces have computed only the leading-order contribution to the quasiparticle interaction 



exactly, while approximately summing certain classes of diagrams to all orders |9l-ll3l|. In par- 
ticular, the summation of particle-particle ladder diagrams in the Brueckner G-matrix was used 
to tame the strong short- distance repulsion present in most realistic nuclear force models, and 
the inclusion of the induced interaction of Babu and Brown [9[ (representing the exchange of 
virtual collective modes between quasiparticles) was found to be essential for achieving the 
stability of nuclear matter against isoscalar density oscillations. 

To date, few works have studied systematically the order-by-order convergence of the quasi- 
particle interaction using realistic models of the nuclear force. In ref. [l^] the pion-exchange 
contribution to the quasiparticle interaction in nuclear matter was obtained at one-loop or- 
der, including also the effects of 2-7r-exchange with intermediate A-isobar states. In the present 
work we derive general expressions for the second-order quasiparticle interaction in terms of the 
partial wave matrix elements of the underlying realistic nucleon-nucleon (NN) potential. The 
numerical accuracy of the second-order calculation in this framework is tested with a scalar- 
isoscalar-exchange potential as well as a (modified) pion-exchange interaction, both of which 
allow for exact analytical solutions at second order. We then study the Idaho N 3 LO chiral 



NN interaction 
interactions 16 
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and derive from this potential a set of low-momentum nucleon-nucleon 
17] . which at a sufficiently coarse resolution scale (A ~ 2fm~ 1 ) provide a 



model-independent two-nucleon interaction and which have better convergence properties when 



employed in many-body perturbation theory [2|, [18| . We extract the four components of the 
isotropic (L = 0) quasiparticle interaction of which two are related to the nuclear matter in- 
compressibility /C and symmetry energy (3. The L = 1 Fermi liquid parameters, associated with 
the angular dependence of the quasiparticle interaction, are used to obtain properties of the 
quasiparticles themselves, such as their effective mass M* and the anomalous orbital ^-factor. 
Our present treatment focuses on the role of two-nucleon interactions. It does not treat the 
contribution of the three-nucleon force to the quasiparticle interaction but sets a reliable frame- 
work for future calculations employing also the leading-order chiral three-nucleon interaction 
In the present work, we therefore seek to identify deficiencies that remain when only 
two-nucleon forces are included in the calculation of the quasiparticle interaction. 

The paper is organized as follows. In Section [Til we describe the microscopic approach 
to Landau's Fermi liquid theory and relate the L = and L = 1 Landau parameters to 
various nuclear matter observables. We then describe in detail our complete calculation of the 
quasiparticle interaction to second order in perturbation theory. In Section II III we first apply 
our scheme to analytically-solvable model interactions (scalar-isoscalar boson exchange and 
modified pion exchange) in order to assess the numerical accuracy. We then employ realistic low- 
momentum nucleon-nucleon interactions and make contact to experimental quantities through 
the Landau parameters. The paper ends with a summary and outlook. 



II. NUCLEAR QUASIPARTICLE INTERACTION 
A. Landau parameters and nuclear observables 

The physics of 'normal' Fermi liquids at low temperatures is governed by the properties and 
interactions of quasiparticles, as emphasized by Landau in the early 1960's. Since quasiparticles 
are well-defined only near the Fermi surface (\p\ = hp) where they are long-lived, Landau's 
theory is valid only for low-energy excitations about the interacting ground state. The quantity 
of primary importance in the theory is the interaction energy between two quasiparticles, which 
can be obtained by functionally differentiating the ground-state energy density twice with 
respect to the quasiparticle densities: 



S 2 S 

J r (piSit 1 ;p 2 s 2 t 2 ) = c~ 
onidn2 



(1) 

ni=ri2=0 



where = ±1/2 and tip = ±1/2 are spin and isospin quantum numbers. The general form of 
the central part of the quasiparticle interaction in nuclear matter excluding tensor components, 



etc., is given by 

Fipufii) = f(puP2) + f(pi,P2)n ■ f 2 + [g(pi,p 2 ) + g'{piiP2)Ti ■ f 2 ] d x ■ B 2 , (2) 

where <ti )2 an d ^1,2 are respectively the spin and isospin operators of the two nucleons on the 
Fermi sphere \p x \ = \p 2 \ = kp. For notational simplicity we have dropped the dependence 
on the quantum numbers Si )2 and t\ >2 , which is introduced through the matrix elements of the 
operators: d\-a 2 — >■ 4sis 2 = ±1 and T\-r 2 — > 4tit 2 = ±1- As it stands in eq. (j2J), the quasiparticle 
interaction is defined for any nuclear density p = 2kp/3iT 2 , but the quantities of physical interest 
result at nuclear matter saturation density p — 0.16 fm -3 (corresponding to kp = 1.33 fm" 1 ). 
For two quasiparticles on the Fermi surface \pi \ — \p 2 \ — kp, the remaining angular dependence 
of their interaction can be expanded in Legendre polynomials of cos# = pi ■ j3 2 : 

X(p 1 ,p 2 ) = J2 X lPl(cos6), (3) 

L=0 

where X represents /, /', g, or g', and the angle 9 is related to the relative momentum p = 
\ Wi ~ P2 1 through the relation 

Q 

p = k F sin - . (4) 

It is conventional to factor out from the quasiparticle interaction the density of states per unit 
energy and volume at the Fermi surface, iVo = 2M*kp/7r 2 , where M* is the nucleon effective 
mass (see eq. ([6])) and kp = 1.33 fm" 1 . This enables one to introduce an equivalent set of 

dimensionless Fermi liquid parameters Fl,Gl, F' l , and G' L through the relation 

00 

= tt J2 [Fl + F' L r x ■ r 2 + (G L + G' L n ■ t 2 )B x ■ a 2 \ P L (cos6). (5) 

Provided the above series converges quickly in L, the interaction between two quasiparticles on 
the Fermi surface is governed by just a few constants which can be directly related to a number 
of observable quantities as we now discuss. 

The quasiparticle effective mass M* is related to the slope of the single-particle potential 
at the Fermi surface and can be obtained from the Landau parameter f\ by invoking Galilean 
invariance. The relation is found to be 

M* M N 3tt 2/1 ' 1 } 

where = 939 MeV is the free nucleon mass. The compression modulus JC of symmetric 
nuclear matter can be obtained from the isotropic (L = 0) spin- and isospin-independent 
component of the quasiparticle interaction 

k=|!(i+fo). m 



The compression modulus of infinite nuclear matter cannot be measured directly, but its value 
/C = 250 ± 50 MeV can be estimated from theoretical predictions of giant monopole resonance 



energies in heavy nuclei 2014221] . The nuclear symmetry energy can be computed from the 



isotropic spin-independent part of the isovector interaction: 

Global fits of nuclear masses with semi-empirical binding energy formulas provide an average 
value for the symmetry energy of = 33 ± 3 MeV over densities in the vicinity of saturated 



nuclear matter 
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24|. The quasiparticle interaction provides also a link to the properties of 
single-particle and collective excitations. In particular, the orbital g-factor for valence nucleons 
(i.e., quasiparticles on the Fermi surface) is different by the amount 5gi from that of a free 
nucleon: [7|: 

1 + r 3 F[ - F 1 _ 1 + r 3 
91=—+ 6(1 + Fl/3 ) r3 = — + ( 9 ) 

One possible mechanism for the anomalous orbital g-factor Sgi are meson exchange currents 



25 



261 ]. which arise in the isospin-dependent components of the nucleon- nucleon interaction. 



According to eq. (Q, the renormalized isoscalar and isovector orbital g-factors are 



(is) 1 / (p) ( n )\ 1 



are different. The former receives no correction, while the latter is sizably enhanced by the 
(reduced) effective mass M* as well as by the (positive) Landau parameter F[. It receives a 
large contribution from one-pion exchange. 

Nuclear matter allows for a rich variety of collective states, including density (breathing 
mode), spin (magnetic dipole mode), isospin (giant dipole mode), and spin-isospin (giant 
Gamow- Teller mode) excitations. As previously discussed, the breathing mode is governed 
by the incompressibility /C of nuclear matter |2jJ. The energy of the (isovector) giant dipole 
mode is correlated with the nuclear symmetry energy f3 j^], while the dipole sum rule 26] 
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dua(El) = ^^(l + 2S 9l ) (11) 



is connected to the anomalous orbital g-factor 5gi with a = 1/137. Experimental results 
are consistent with a value of the anomalous orbital g-factor of 5gi ~ 0.23 ± 0.03. Finally, 
the giant Gamow- Teller resonance has been widely studied due to its connection to the nuclear 
spin-isospin response function and for ruling out pion condensation in moderately-dense nuclear 



matter. An analysis of the experimental excitation energies and transition strengths 291431] 
leads to a value for the parameter 

0^-0.6-0.7, (12) 
which is used to model the spin-isospin interaction in nuclei as a zero-range contact interaction. 
As a convention it is related to the dimensionless Landau parameter G' by 

9 2 

Go = No ^2 9nn j ( 13 ) 
where — 13.2 is the strong tcN coupling constant. It is well-known that the giant Gamow- 



Teller resonances receive important contributions from the coupling to A-hole excitations 31]. 



Such dynamical effects due to non-nucleonic degrees of freedom are reflected in the leading- 



order, 27r-exchange three-nucleon interaction to be included in future work 19]. 



B. Quasiparticle interaction at second order 

Expanding the energy density to second-order in the (Goldstone) linked-diagram expansion 
and differentiating twice with respect to the nucleon distribution function, one obtains for the 
first two contributions to the quasiparticle interaction 

T {1 \p lSl tup 2 s 2 t 2 ) = (pis 1 t 1 ;p 2 s 2 t 2 \V\p 1 s 1 t 1 ;P2S2t 2 ) = (12|F|12) (14) 

and 

|<12|F|34)| 2 (l-n 3 )(l-n 4 ) 



J r(2 \piS 1 t 1 ;p2S2t 2 ) = -J 



2 ^ ei + e 2 - e 3 - e 4 



|<12|t/|34}| 2 n 3 n 4 | (13|K|24) | 2 n 3 (l - n 4 



2 ^ e 3 + e 4 - ex - e 2 ^ e x + e 3 - e 2 - e 4 
In eqs. (|14|) and (|15p the quantity V denotes the antisymmetrized two-body potential (with units 
of fm 2 ) given by (plSJT\V\p'l'SJT) = (1 - (-l) L+s+T )(plSJT\V\p'l'SJT) in the partial wave 
basis, and in eq. (ITS]) the summation is over intermediate-state momenta, spins and isospins. We 
specify our sign and normalization conventions through the perturbative relation between di- 
agonal two-body matrix elements and phase shifts: tanSisj(p) = —pM^{plS 'JT\V\plS 'JT) j 'Am . 
The first-order term of eq. (TT4l is just the diagonal matrix element of the antisymmetrized 
two-body interaction, while the second-order term (eq. (j!5p ) has been separated into particle- 
particle, hole-hole, and particle-hole terms depicted diagrammatically in Fig. [TJ The distribu- 
tion function n is the usual step function for the nuclear matter ground state: 

, 1 for \k\ < kp 

W~ ■ ( l6 ) 

for \k | > kp 
6 






FIG. 1. Diagrams contributing to the second-order quasiparticle interaction (exchange terms omitted): 
(a) particle-particle diagram, (b) hole-hole diagram, and (c) particle-hole diagram. 

In the following, we discuss the general evaluation of eqs. (|14p and f[To"j) for interactions given 
in the partial- wave basis. We first define the spin-averaged quasiparticle interaction F(pST), 
which is obtained from the full quasiparticle interaction by averaging over the spin-substates: 



(17) 



where p = — p<± \ = kpsm.6/2, and in ^{pSmsTT;,) the spins and isospins of the two 
quasiparticles are coupled to total spin S = 0, 1 and total isospin T = 0, 1. We take an isospin- 
symmetric two-body potential and thus the quasiparticle interaction is independent of T z . The 
first-order contribution to the central part of the quasiparticle interaction is then obtained by 
summing over the allowed partial wave matrix elements: 



? {1 \pST) = Y1(2J + 1) (plSJT\ V \plSJT) . 



;is) 



Note that there is an additional factor of 4tt in eq. (41) in ref. 32J and eq. (28) in ref. 13| due to a 
different normalization convention. From eq. (fl~8]) we can project out the individual components 
of the quasiparticle interaction using the appropriate linear combinations of F(pST) with S = 
0,1 and T = 0,1: 

f(p) = ^ (1) (p00) + ^ (1) (p0i) + ^ {1 \ P io) + ^T (1) ( P n) 
g(p) = -^ (1) (p00) - ^ (1) (pQi) + ^ {1 \ P io) + ^ (1) (pn) 
fip) = -^ (1) (p00) + ^ (1) (p0i) - ^ (1) (pio) + ^ (1) (pn) 
g'ip) = ^ {1 \p00) - y/ {1 \p^) - ^ (1) (pio) + ^ (1) (pH). (19) 

The leading-order expressions, eqs. ({TBI) and (|T9l . give the full p-dependence (i.e., angular 
dependence) of the quasiparticle interaction, and therefore one can project out the density- 
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dependent Landau parameters for arbitrary L: 

X L = 2(2L + 1) f F dp^X{p)P L {l - 2p 2 /k 2 F ). (20) 

JO fc F 

For the second-order contributions to the quasiparticle interaction, the complete p-dependence 
is in general not easily obtained (e.g., for the particle-hole term). We instead compute the 
Landau parameters for each L separately, choosing the total momentum vector to be aligned 
with the z-axis. In the following, the two quasiparticle momenta are labeled p\ and p 2 , while 
the intermediate-state momenta are labeled k 3 and k\. For the particle-particle contribution 
one finds 

jf pp (Sm s T) = Yl J dpP J dqq 2 N(l im l 2 m%mkm>)P™(0)P™(0) 

m'JJ'M 



Mjq i i rmm{x ,li f ^ 

X o _ n 2 * 2 3 1 4 ^hmsm a ^hrn'sm' s ^l3msm a ^Um'sm' s / C? COS 9 q P^ [COS 9q)P™ [COS 9q) 

P y J max{— xo,— 1} 

x{pl 1 SJMT\V\ql 2 SJMT){qUSJ'MT\V\pl z SJ'MT)P L {l - 2p 2 /k 2 F ), (21) 



where P™ are associated Legendre functions, p — (pi — P2)/2, Q — (&3 ~~ ^0/2, xq = (q 2 



p 2 )/2q^k 2 F - p 2 and N (hml 2 m' hml 4 m') = N%N%'N™N™' with Nf 1 = y/(2l + I) {I - m)!/(/ + m) 
Similarly, for the hole-hole diagram one obtains 

Tf hh (Sm s T) = Y, d pp dqq 2 N(hmkm%mhm')P™(0)PZ(0) 

F 7 7 7 7 / ** */ 

m' s JJ'M 

M N i i i i JM j M ji M i /•mia{-x ,l} 

X 2 _ 2^ 2 3 1 4 ^hrnsm a ^l2m'sm l J^l 3 msm a ^l i rn'sm' a I d COS 9qPi 2 (COS 9 q )Pi 4 (COS #<j) 

9 P Jmax{xo,- 1} 

x(pl 1 SJMT\V\ql 2 SJMT)(qUSJ'MT\V\phSJ'MT)P L {l - 2p 2 /k 2 F ). (22) 

Averaging over the spin substates and employing eq. ( fT9l) with the substitution J 7 *- 1 - 1 — >• J^ 2 -* 
again yields the individual spin and isospin components of the quasiparticle interaction. The 
evaluation of the particle-hole diagram proceeds similarly; however, in this case the coupling 
of the two quasiparticles to total spin (and isospin) requires an additional step. Coupling to 
states with m s = is achieved by taking the combinations (neglecting isospin for simplicity) 

^ (2)p ViP 2 ; S = 1/2 ± 1/2, = 0) = -2^2 [(Pih t s 3 \V\p 2 k 4 I s 4 )(p 2 h I s 4 |V|&£ 3 t s 3 ) 

34 

± (Pih t s 3 \V\p 2 k 4 t s 4 )(p 2 k 4 I s 4 |VK£ 3 I s 3 )} ^ (1 ~ n4) ■ (23) 



ei + e 3 - e 2 - e 4 

We provide the expression for J 7 ^ 2 ^ (corresponding to L = 0) in uncoupled quasiparticle spin 
and isospin states, appopriate for evaluating the first term in eq. (123]) . which can be easily 



generalized in order to obtain the second term. We find 



I r-2k F ry k F -P / 4 rk F +p'/2 

J r Q 2 ^ ph {siS2tit 2 ) = / dp' / dpp / dqqN^imhm'lzmUm') 

^ k F .stfafc J ° Jk F - P >/2 J^W^l 

lil^lsl&mm' 
SJMTS'J'T' 

_^N_pm( n )P m (m^ff )P m '(m*fl )P m '(rn^f) \jh+h-h-hnJM nJM nJ'M nJ'M 

xC(s 1 s 2 S3S4)C(t 1 t 2 t3t4)(^i5JMT|V|g/ 2 SJMr)(g/45V'Mr , |V|p/3S'J , MT'), (24) 

where now p = (pi — fca)/2, q = (p 2 — &4)/2, and the total momentum p ' = p\ + = p 2 + 
The angle 8 P between p and p' is fixed (via \p\ \ = kp together with p\ = p' /2 + p) by the 
relation pp' cos 6 P = k 2 F — p 2 — p' 2 /4, and analogously for the angle 9 q between q and p'. The 
combination of spin Clebsch-Gordan coefficients that arises in the above expression is denoted 
by 

C(s lS2 s 3 s 4 ) = Cf m \ Cf m \ Cf m r Cf m i , (25) 

V ' 2 S l2 S 3 2 S 2^4 2 S l2 S 3 2 S 22 S 4 

and likewise for the combination of isospin Clebsch-Gordan coefficients. In computing the 
particle-hole term for L > 0, we use 



Pl{p\ - fa) = Pl (jp-P' 2 + 2p"P'(pcos0j> + gcos0 9 ) + || cos O > 



(26) 



and employ the addition theorem for spherical harmonics to write cos 6^ in terms of cos8 p , 
cos q and an azimuthal angle <fi. The involved integral d<j) gives different selection rules for 
the m,m' values of the associated Legendre functions in eq. fl24l . In deriving eqs. fl2T]) - fl24|) . 
we have assumed that the intermediate-state energies in eq. (JIB"]) are those of free particles: 
£/c = k 2 /2M]y. Later we will include the first-order correction to the dispersion relation arising 
from the in-medium self-energy, which leads to the substitution — > M* in the above 
equations. 



III. CALCULATIONS AND RESULTS 



A. Prelude: One boson exchange interactions as test cases 

The numerical computation of the quasiparticle interaction at second order is obviously quite 
intricate, and truncations in the number of included partial waves and in the momentum- space 
integrations are necessary. In such a situation it is very helpful to have available analytical 
results for simple model interactions in order to test the accuracy of the numerical calculations. 
For that purpose we derived in this subsection analytical expressions for the quasiparticle 
interaction up to second order arising from (i) massive scalar-isoscalar boson exchange and 
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(a) (b) 



(c) 



(d) (f) 



(g) 



FIG. 2. Diagrammatic contributions to the second-order quasiparticle interaction in nuclear matter 
organized in the number of medium insertions (symbolized by the short double lines). Crossed dia- 
grams and additional reflected diagrams are not shown. The labels (b) and (f) refer to the crossed 
terms of (a) and (d). 

(ii) pion exchange modified by squaring the static propagator. We omit all technical details 



of these calculations which can be found (for L = 0) in ref. 14j for the case of tree-level 
and one-loop (i.e., second-order) pion-exchange. In the present treatment the second-order 
quasiparticle interaction is organized differently than in eq. ( 115)) . The explicit decomposition of 
the in-medium nucleon propagator into a particle and hole propagator is replaced by the sum 
of a "vacuum" and "medium insertion" component: 

9(\p\-k F ) 0(k F -\p\) 



G(po,p) 



p - p 2 /2M N + te + po- P 2 /2M N - te 



-2n5(p -p 2 /2M N )9(k F -\p\), 



-..) ,-, i r - -'•'•>/•>» /• / \ /■■ \' - r ir i '• (27) 

Po - p 2 /2M N + te 

and the organization is now in the number of medium insertions rather than in terms of particle 
and hole intermediate states. The central parts of the quasiparticle interaction are constructed 
for any L through an angle-averaging procedure 

2L + 1 



(Any 



dQ.xdQ 2 {piP2\V&\piP2)PL{vi - fa) 



/l(^f) + 9L{kF)vi ■ <?2 + f' L (kF)ri ■ r 2 + g' L (k F )at ■ a 2 n ■ f 2 



(2f 



In this equation V e g represents the effective model interaction computed up to one-loop order 
(second order). 

We first consider as a generic example the exchange of a scalar-isoscalar boson with mass 
m s and coupling constant g s (to the nucleon). In momentum and coordinate space it gives rise 
to central potentials of the form 



V C {q) 



91 



mj + q 2 



V c {r) 



9 2 s e- m ° r 
Att r 



For the first-order contributions to the L = 0, 1 Landau parameters one finds 



4 1} (k F ) = ^ 



9; 



mi 



1 + (1 + <t)(1 + t 



ln(l + 4u 2 
16u 2 



(29) 



(30) 
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;i + «r)(l + T) 



32m 2 w 4 



l + 2w 2 )ln(l + 4w 2 



4m z 



(31) 



where cr = d\ ■ a% and r = f\ ■ T2 are short-hand notations for the spin-spin and isospin-isospin 
operators. The dimensionless variable u = kp/m s denotes the ratio of the Fermi momentum kp 
to the scalar boson mass m s . Note that in this approach both direct and crossed diagrams can 
contribute. The crossed diagrams have to be multiplied by the negative product of the spin- 
and isospin-exchange operators —(1 + c)(l + r)/4. At second order there are five classes of 
diagrammatic contributions, shown in Fig. [2} to the quasiparticle interaction. The direct terms 
from iterated (second order) boson exchange, see Fig. Eta), read 



4 2a) (k F ) 



g*M N ln(l + 4w 2 ) 



F? a \k F ) 



3g*M N 



327rm^ 

Au 2 - (1 + 2u 2 ) ln(l + 4m 2 ) 



(32) 
(33) 

(34) 

167rm3 J u\l + 2x 2 ) L — ~j • ( 35 ) 
The coupling of the exchanged boson to nucleon-hole states, Fig.^c), gives rise to nonvanishing 
crossed terms which read 



647rm^M 4 

whereas the corresponding crossed terms (b) have the form 

arctan 2x — arctan x 



4 2b \kp) = (l+CT)(l 

4 2b \k F ) = (l + a)(l + T 



gjM N 

167rmjj jq 



dx 



3^M 



N 



dx 



u 



2x 2 



u 2 (l + 2x 2 ) 

■ [ arctan 2x — arctan x] 



(2c), 



(k F ) = (1 + + t) 



N 



An 2 m 3 u jq 



dx ■ 



1 



(1 + <t)(1 + t) 



3g*M N 



dx 



(l + 4x 2 ) 2 
u 2 - 2x 2 



2ux + (u 2 — x 2 ) In 
2ux + (u 2 — x 2 ) In 



u + x 
u — x 

u + x 
u — x 



(36) 



(37) 



A7r 2 m 3 s u 4 J (1 + Ax 2 ) 2 
Pauli blocking occurs in the planar- and crossed-box diagrams, Fig. |2fd)-(e), and for the sum 

of their direct terms one finds the forms 







• ^(kp 
(k F ) 



N 



ir 2 m 3 u 2 jq 



dx ■ 



(2d+2e) . 



gjM N 

7i 2 m 3 u 4 jq 



dx 



1 + Ax 2 ) 2 



x 

l + 4x 2 ) 2 



u In • 



u 



u — X 



x In 



u — x 

x 



(3f 



2ux(x — u) + 2m 3 In ■ 



u 



+3x(u 2 — x 2 ) In 



u 



X 



+ 2x A In ■ 



u — X 
u — X 



(39) 



On the other hand, the crossed terms of the planar-box diagram with Pauli blocking, see Fig. 
W) yield 

.(2/), q 4 *Mn r. 



n n {kp) 



l + -)(l + r) 



4ii 2 m 3 u 2 Jq 



dx 



l + Ax 2 



"* J_ ln uy/R+(l-4xy)(x-y) 
VR uVR+ (Axy - l)(x - y) 



(40) 
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;i+<r)(i+r) 



N 



. . . 2 1 + {U + Xf 

dx{ W ) '— + 



u In 



4(u — x) 



-x In 



u 



x 



dy 



l + 2u 2 



l + (u-x) 2 1 + Ax 2 

1 — 4xy)(x — y) 



u + x 



Ax 2 , mV^R 
In 



(41) 



'R uVR+{Axy-l)(x-y)_ 

with the auxiliary polynomial R = 4u 2 + (4x 2 — 1)(% 2 — 1). Finally, the density-dependent 
vertex correction to one-boson exchange, Fig. [2(g), provides a nonzero contribution only in the 
crossed diagram. The corresponding expressions for the L — 0,1 Landau parameters read 

ln(l + Ax 2 



;i+<r)(i- 



16n 2 m^u 



N 

3/1,2 



dx 



\/l + 4m 2 - Ax 2 



v u\/l + 4w 2 - 4x 2 + x) 2 

in — ; Trz~, - ^ , (,4^J 



l + 4M 2 )(M 2 -a; 2 ) 



jf\k F ) = (1 + <r)(l + r)-|4^fl f^ 1 ^ 1 + ln — 

o2tT TTlgU Jq { u — X 

l + 2u 2 , (u v / 1 + 4m 2 -4x 2 + x) 2 1 

ln ^ i TT^T^ „^ r • ( 43 ) 



VI + 4m 2 - 4x 2 (1 + 4m 2 )(m 2 -x 2 

Since at most double integrals over well-behaved functions are involved in the expressions in 
eqs. (15U|) - (I43I) . they can be evaluated easily to high numerical precision. After summing them 
together, they provide a crucial check for our calculation of the second-order quasiparticle 
interaction in the partial wave basis (see Section III Bp . We set the scalar boson mass m s = 500 



MeV and coupling constant g s = 2.5 and work with the partial wave matrix elements following 
from the central potential Vc(q) in eq. (l29l) . 

Table [J shows the dimensionful Fermi liquid parameters (labeled 'Exact') as obtained from 
the above analytical formulas at nuclear matter saturation density (kp = 1.33 fm _1 ). Due to the 
simple spin and isospin dependence of the underlying interaction the constraint 9l = fi = 9l 
holds. For comparison we show also the first- and second-order results obtained with the general 
partial wave expansion. The second-order terms are further subdivided into particle-particle, 
hole-hole and particle-hole contributions. We find agreement between both methods to within 
1% or better for all L = 0, 1 Fermi liquid parameters. In order to achieve this accuracy, the 
expansions must be carried out through at least the lowest 15 partial waves. 

A feature of all realistic NN interactions is the presence of a strong tensor force, which results 
in mixing matrix elements between spin-triplet states differing by two units of orbital angular 
momentum. At second order these mixing matrix elements generate substantial contributions 
to the L = 0, 1 Fermi liquid parameters. In order to test the numerical accuracy of our partial 
wave expansion scheme for the additional complexity arising from tensor forces, we consider 
now the quasiparticle interaction in nuclear matter generated by (modified) "pion" exchange. 
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Scalar-isoscalar boson exchange {kp = 1.33 fm 1 ) 




fo [fm 2 ] 


9o [fm 2 ] 


f'o If™ 2 ] 


g'o [fm 2 ] 


fx [fm 2 ] 


9i [fm 2 ] 


f> [fm 2 ] 


g'i [fm 2 ] 


1st 


-0.809 


0.164 


0.164 


0.164 


0.060 


0.060 


0.060 


0.060 


2nd(pp) 


-0.186 


0.056 


0.056 


0.056 


0.038 


-0.006 


-0.006 


-0.006 


2nd(hh) 


-0.033 


0.010 


0.010 


0.010 


0.042 


-0.013 


-0.013 


-0.013 


2nd(ph) 


0.198 


0.061 


0.061 


0.061 


0.100 


0.085 


0.085 


0.085 


Total 


-0.830 


0.291 


0.291 


0.291 


0.240 


0.127 


0.127 


0.127 


Exact 


-0.830 


0.292 


0.292 


0.292 


0.242 


0.127 


0.127 


0.127 



TABLE I. Fermi liquid parameters {L = 0, 1) for a scalar-isoscalar boson-exchange interaction with 
parameters described in the text. The exact results at kp = 1.33 fm -1 obtained from our derived 
analytical expressions are compared to the numerical results computed via a partial wave expansion. 



To be specific we take a nucleon-nucleon potential in momentum space of the form 



V T (q) 



r 



[m 2 + q 



^ eri • q a 2 ■ q n • r 2 



(44) 



where g is a dimensionless coupling constant and m n a variable "pion" mass. The isovector 
spin-spin and tensor potentials in coordinate space following from Vr(q) read 



Ve(r) 



g 2 e~ m * r 



(m n r - 2) , 



Vr(r) 



g 2 e- m * r 



1 + m^r) . 



(45) 



24tt r v " w 24vr r 

The basic motivation for squaring the propagator in eq. f )44|) is to tame the tensor potential 
at short distances, and thereby one avoids the linear divergence that would otherwise occur in 
iterated (second-order) one-pion-exchange. In the presence of non-convergent loop integrals, 
analytical and numerical treatments become difficult to match properly. Let us now enumerate 
the contributions at first and second order to the L = 0, 1 Landau parameters as they arise 
from modified "pion" exchange. 
The first-order contributions read 



Jo (1) (M = (3-<r)(3 



J- 1 (1) (M = (3-^)(3-r) ir 



12m 2 L 



1 



4u 2 



4m 2 



u 



ln(l + An 2 



1 



1 



u 



2 l + Au 2 



(46) 
(47) 



with the abbreviation u = kp/m^. For the second-order contributions we follow the labeling 
(a) — (g) introduced previously for scalar-isoscalar boson exchange: 



4 2a \k F ) = (3 - 2t 



N 



3 + 16u 2 
3(1 + Au 2 ) 2 
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ln(l + 4m 2 



(4£ 



F {2a \k F ) = (3-2r) 



9 4 M N 
167rm^. 



35 5 + 24m 2 



8u 2 (l + 4w 2 ) 2 32u 4 



(7 + 6m 2 ) ln(l + 4m 2 



(49) 



4 2b \k F ) = (5r-3) 



g A M N f 1 + cr 



327rm3 \8(1 + m 2 
I F. l + 4x 2 + 8x 4 



3) 



1 , l + 2w 2 

n I 11 7T 

2u 2 1 + it 



1 + 2m 2 



dx 



J**\k F ) = (5t - 3) 



9 4 M N f3 



(<r-3) 



+ dx 



327rm3 1 4 
1 



w 2 (l + 2x 2 ) 3 
(1 + cr 



(arctan2x — arctanx) 



ln(l + u 2 



_l + 2w 2 // 

u 2 - 2x 2 



u 2 2(1 + u 2 ) u 4 

1 2 + w 2 l + 2w 2 
+ — ^-r- In 



u 4 (l + 2x 2 ) : 



2 2m 4 1 + u 2 
(1 + 4x 2 + 8x 4 )(arctan2x — arctanx) 



^ 2c) (M = (3-^)(3-r) 



% 4 M, 



N 



dx 



x 



3 r K 2 m\u 2 J (1 + 4x 2 ) 



n 2 rri^u 4 

-A 



T (2c) ( , v w, ^ r X> 2 -2x 2 ) 

5 



lux + (w 2 — x 2 ) In - X 
it — X 



2ux + (u — x ) In 

it — x 



^ ) (M = (3-2r) 16 ^ Mj ! / dx ' 



7r 2 m^ 2 7 (l + 4x 2 ) z 



U + X u 2 — x 2 
u In — r + x In — 



4(w — x) 



x z 



(50) 



(51) 

(52) 
(53) 
(54) 



^ 2d \k F ) = (3-2t) 



l(i//'.\/v / , 



dx 



7r 2 m3 M 4 7 (l + 4x 2 ) 



2m" 1 (w - x) + W±n 



31. u + x 



4(u — x) 



+x(3u 2 - 2x 2 ) In 



u 2 — x 21 



^ 2e) (M = (3 + 2r) 



16ff 4 Mjv r, a; 5 



dx ■ 



Ti 2 mtu 2 7n (1 + 4x 2 ) 4 [ u 2 — x 2 



Au 2 

u In — — x In 



Ff e \k F ) = (3 + 2r) 



16(/ 4 Mjv F 



dx 



x 



TT 2 mlu 4 ./n (l + 4x 2 ) 4 



2u(x 2 -u 2 )+m 3 In 



4m 2 



■u 2 — X 2 



x z 

u + x 
u — X 

s , M + X 

•x In 



u — X 



(55) 

(56) 
(57) 



We split the crossed terms from the planar-box diagram with Pauli blocking, see Fig. EJ^f), into 
factorizable parts: 

2 



F (2f \k F ) = (3+<t)(3-5t) _ g * M " „ / r/.< 



J-f^(A; F ) = (3 + o-)(3-5r) 



967r 2 m 3 M 2 J 



x — u 



X + u 1 , 1 + (m + x) 2 
H In 



1 + (u — x) 2 1 + (u + x) 2 2x 1 + (u — x) 



g 4 M N 



9Qn 2 m^u 2 



/ dx 


G 2 + 2G 2 " 


/o 





M + X X — U 

(^a — ~ — : — - — : — — + 



1 1 + (u + x) 2 
1 + (m + x) 2 ' l + {u-x) 2 2u n i + (u-x) 2 ' 

3 ^ 3 /, 2 2m 1 + (U + X) 2 



(58) 
(59) 

(60) 

(61) 
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and non-factorizable parts: 



T 



(2/0, 



2g A M h 



<i • - - ^ r;, x " r ",, f M^-y)(i+4^) 

„ (M = (--3)(3-5r ) ^^y o ^ ITTi - ¥ y o ^( fl[4(a; _ y)2 + i2] 



1/9 9 9 „ 9 9n i mVR + (1 — Axy)(x — y) 
(m 2 - x 2 - + 8a; V) In 



W>l 2 



m 



Vft + (4xy- l)(x-y)J ' 



(62) 



jf f '\k F ) = (<r-3)(3-5r) 



647r 2 m 3 r M 4 jg 



. 1 + (m + x) 2 
In- , + 



1 + (u - x) 2 



;i + 4x 2 ) 2 



m In 



1 + (it - x) 2 1 + (it + x) 2 
4(m — x) 



u + x 



x In 



m 2 — x 2 



+ 



'Au(x-y)(l + Axy)(l + 2u 2 -Ax 2 ) 

i)[4(,-y) 2 + ^] ^ +f2(1 H2 "^ 



x ^ _ 2 _ 2 8 V) fl \ J_ ln uVR+{l-Axy){x-y) \ 1 1 

with auxiliary polynomial i? = 4m 2 + (4x 2 — l)(4y 2 — 1). These two pieces, (2/) and (2/'), are 
distinguished by whether the remaining nucleon propagator can be cancelled or not by terms 
from the product of (momentum-dependent) ttN interaction vertices in the numerator. Finally, 
the density- dependent vertex corrections to modified "pion" exchange have nonzero crossed 
terms, which we split again into factorizable parts: 

Au 2 



^ 9 \k F ) = {3 - a){3 - r) 



9 A M N 
967i 2 miu 3 



1 + Au 2 



ln(l + 4m 2 ) 



l-i±^ln(l + 4 M 2 ) 
4m 2 



(64) 



F? 9 \k F ) = {3 - *){3 - r) 



<? 4 M 



N 



327r 2 m 3 r M 5 



1 - 



1 + 2m 2 
4m 2 



ln(l + 4m 2 ) 



x 



3m 2 + 



u 



1 + 4m 2 



(1 + M 2 )ln(l + Au 2 ) 



and non-factorizable parts: 

4 29 '\k F ) = (3 - <r)(3 - t) 
2ux(l + Au 



N 



2Aix 2 m\u 2 jo 



dx 



Ax' 



l + Ax 2 



-ln(l + 4x 2 



2\-l 



V 2 — X 2 



1 + 4m 2 -4x 2 (1 + 4m 2 - 4x 2 ) 3 / 2 



In 



u 



VI + 4m 2 - 4x 2 



x) 



'I + Au 2 ){u 2 - x 2 ) 



(65) 



(66) 



F? 9 '\k F ) = (3 - <x)(3 - r 



9 4 M N 



x 



327r 2 m> 4 jo 
Aux(l + 2u 2 ) 
(1 + Au 2 ){l + Au 2 - Ax 2 ) 



dx 



Ax z 



-ln 



1 + 4x 2 

M + X 



ln(l + Ax 2 ) 



u — X 



x 2 ) (6 + 4m 2 ) (wVl + 4m 2 - Ax 2 + x) 2 



(1 + 4m 2 -4x 2 ) 3 / 2 



In 



(1 + 4m 2 ) (m 2 - x 2 ) 



(67) 
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Modified "pion" exchange (kp = 1.33 fm 1 ) 




fo [fm 2 ] 


9o [fm 2 ] 


f'o If™ 2 ] 


g'o [fm 2 ] 


fx [fm 2 ] 


9i [fm 2 ] 


f> [fm 2 ] 


9'i [fm 2 ] 


1st 


0.244 


-0.081 


-0.081 


0.027 


-0.079 


0.026 


0.026 


-0.009 


2nd(pp) 


-0.357 


-0.062 


0.269 


0.104 


0.018 


-0.005 


0.027 


0.009 


2nd(hh) 


-0.017 


-0.002 


0.009 


0.003 


0.029 


0.003 


-0.014 


-0.005 


2nd(ph) 


0.146 


-0.023 


0.027 


0.008 


0.008 


0.010 


0.036 


-0.003 


Total 


0.017 


-0.169 


0.224 


0.142 


-0.024 


0.035 


0.075 


-0.009 


Exact 


0.017 


-0.169 


0.224 


0.142 


-0.023 


0.035 


0.074 


-0.009 



TABLE II. Fermi liquid parameters (L = 0, 1) for (modified) pion exchange. The exact results at 
kp = 1.33 fm -1 are obtained from our derived analytical expressions and compared to the numerical 
results computed via a partial wave expansion. 

Together with the coupling constant g = 2.5 we choose a large "pion" mass m n = 400 MeV in 
order to suppress partial wave matrix elements from the model interaction Vr(q) beyond J = 6 
in the numerical computations based on the partial wave expansion scheme. 

We show in Table [III the L — 0,1 Fermi liquid parameters (at kp = 1.33 fm -1 ) for the 
modified "pion" exchange interaction up to second order in perturbation theory. The summed 
results from the analytic formulas eqs. fj46l) - fl67j) are labeled "Exact" and compared to the 
results obtained by first evaluating the interaction in the partial wave basis and then using eqs. 
()2ip - (|24p . As in the case of scalar-isoscalar exchange, we find excellent agreement between the 
two (equivalent) methods. 



B. Realistic nuclear two-body potentials 

After having verified the numerical accuracy of our partial wave expansion scheme, we 
extend in this section the discussion to realistic nuclear two-body potentials. We start wit 
the Idaho N 3 LO chiral NN interaction [15j and employ renormalization group methods 



16- 



181 ] to evolve this (bare) interaction down to a resolution scale (A ~ 2 fm 1 ) at which the 



NN interaction becomes universal. The quasiparticle interaction in nuc 



13 



32 



ear matter has been 



331 ] . but a complete 



studied previously with such low-momentum nuclear interactions 
second-order calculation has never been performed. Given the observed better convergence 
properties of low-momentum interactions in nuclear many-body calculations, we wish to study 
here systematically the order-by-order convergence of the quasiparticle interaction derived from 
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low-momentum NN potentials. A complete treatment of low-momentum nuclear forces requires 
the consistent evolution of two- and three-body forces together. We postpone the inclusion of 



contri 



Dutions to the quasiparticle interaction from the (chiral) three-nucleon force to upcoming 



work [19 ]. 



In Table II III we compare the L — 0, 1 Fermi liquid parameters obtained from the bare 
chiral N 3 LO potential to those of low-momentum interactions obtained by integrating out 
momenta above a resolution scale of A = 2.1 fm -1 and A = 2.3 fm -1 . The intermediate- 
state energies in the second-order diagrams are those of free nucleons = k 2 /2M^, and 
we include partial waves up to J = 6 which result in well-converged L = 0,1 Fermi liquid 
parameters. Comparing the results at first-order, we find a large decrease in the isotropic spin- 
and isospin-independent Landau parameter fo as the decimation scale decreases. This enhances 
the (apparent) instability of nuclear matter against isoscalar density oscillations. The effect 
results largely from integrating out some short-distance repulsion in the bare N 3 LO interaction. 
A repulsive contact interaction Vc = 4C (contributing with equal strength 4C in singlet and 
triplet 5-waves) gives rise to a first-order quasiparticle interaction of the form 

Tip = C(3 — &x ■ a 2 — Ti • f 2 - a\ ■ a 2 Ti • r 2 ) (68) 

and no contributions for L > 1. Thus, integrating out the short-distance repulsion in the chiral 
N 3 LO potential yields a large decrease in fo and a (three-times) weaker increase in go, fo, and 
g'o- The increase in /q gives rise to an increase in the nuclear symmetry energy at saturation 
density by approximately 20% for interactions evolved down to A ~ 2 fm -1 . Overall, the scale 
dependence of the first-order L = 1 Landau parameters is weaker, and in particular the two 
isospin-independent (fi and gi) components of the L = 1 quasiparticle interaction are almost 
scale independent. However, the parameter f[ increases as the cutoff scale is lowered, which 
results according to eq. ([9]) in an increase in the anomalous orbital g-factor by 10-15%. 

Considering the three parts comprising the second-order quasiparticle interaction, we find 
large contributions from both the particle-particle {pp) and particle-hole (ph) diagrams. In 
particular, the ph term is quite large, which suggests the need for an exact treatment of this 
contribution which until now has been absent in the literature. As the decimation scale is 
lowered, the pp contribution is generally reduced while the hole-hole (hh) and ph contributions 
are both increased. In previous studies, the hh diagram has often been neglected since it was 
assumed to give a relatively small contribution to the quasiparticle interaction. However, we 
learn from our exact calculation that its effects are non-negligible for all of the spin-independent 
Landau parameters. 
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Idaho N 3 LO potential for k F = 1.33 fin" 1 




fo [fm 2 ] 


90 [fm 2 ] 


fo [fm 2 ] 


9'o [fm 2 ] 


fi [fm 2 ] 


91 [fm 2 ] 


fi [fm 2 ] 


g'i [fm 2 ] 


1st 


-1.274 


0.298 


0.200 


0.955 


-1.018 


0.529 


0.230 


0.090 


2nd(pp) 


-1.461 


0.023 


0.686 


0.255 


0.041 


-0.059 


0.334 


0.254 


2nd(hh) 


-0.271 


0.018 


0.120 


0.041 


0.276 


0.041 


-0.144 


-0.009 


2nd(ph) 


1.642 


-0.057 


0.429 


0.162 


0.889 


-0.143 


0.130 


0.142 


Total 


-1.364 


0.281 


1.436 


1.413 


0.188 


0.367 


0.550 


0.477 


Mow-k(A = 2.3 fin" 1 ) for k F = 1.33 fin" 1 




fo [fm 2 ] 


90 [fm 2 ] 


fo [fm 2 ] 


g'o [fm 2 ] 


fi [fm 2 ] 


91 [fm 2 ] 


fi [fm 2 ] 


g[ [fm 2 ] 


1st 


-1.793 


0.357 


0.394 


1.069 


-0.996 


0.493 


0.357 


0.152 


2nd(pp) 


-0.974 


-0.098 


0.594 


0.185 


-0.129 


-0.003 


0.252 


0.193 


2nd(hh) 


-0.358 


0.030 


0.169 


0.075 


0.338 


0.028 


-0.180 


-0.042 


2nd(ph) 


2.102 


0.095 


0.588 


0.254 


1.512 


0.003 


0.204 


0.329 


Total 


-1.023 


0.385 


1.744 


1.583 


0.725 


0.521 


0.634 


0.632 


Mow-k(A = 2.1 fm- 1 ) for k F = 1.33 fin" 1 




fo [fm 2 ] 


90 [fm 2 ] 


fo [fm 2 ] 


g'o [fm 2 ] 


fi [fm 2 ] 


9i [fm 2 ] 


fi [fm 2 ] 


g[ [fm 2 ] 


1st 


-1.919 


0.327 


0.497 


1.099 


-1.034 


0.475 


0.409 


0.178 


2nd(pp) 


-0.864 


-0.079 


0.507 


0.164 


-0.130 


0.011 


0.236 


0.174 


2nd(hh) 


-0.386 


0.022 


0.195 


0.085 


0.355 


0.034 


-0.195 


-0.049 


2nd(ph) 


2.033 


0.164 


0.493 


0.292 


1.620 


0.098 


0.234 


0.412 


Total 


-1.135 


0.434 


1.692 


1.640 


0.812 


0.617 


0.684 


0.715 



TABLE III. Fermi liquid parameters (L = 0, 1) for the Idaho N 3 LO chiral potential as well as for the 
low-momentum nucleon-nucleon interaction Vj ow -k at a cutoff scale of A = 2.1 and 2.3 fm -1 . 

At second order, the contributions to f are sizable and approximately cancel each other for 
the bare N 3 LO chiral NN interaction, but they become more strongly repulsive as the resolution 
scale A is decreased. This reduces the large decrease at leading-order in / effected through 
the renormalization group decimation, so that after including the second-order corrections, the 
spread in the values of f for all three potentials (bare N 3 LO and its decimations to A = 2.1 fm -1 
and A = 2.3 fm" 1 ) is much smaller than at first order. For each of the three different potentials, 
the second-order terms are strongly coherent in both the f and g' channels. In the former 
case, this change alone would give rise to a dramatic increase the nuclear symmetry energy (3. 
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FIG. 3. Density dependence of the Fermi liquid parameters for the low-momentum NN interaction 
Mow-k(iY = 2.1 fm -1 ). Here po = 0.16 fm -3 is the nuclear matter saturation density. 



This effect will be partly reduced through the increase in the quasiparticle effective mass M*, 
which we see is now close to unity for the bare N 3 LO chiral interaction but which is strongly 
scale-dependent and enhanced above the free mass as the decimation scale is lowered. 
The parameter g' Q , related to the energy of Giant Gamow- Teller resonances, is increased by 
approximately 50% after inclusion of the second-order diagrams. In Fig. [3] we plot the Fermi 
liquid parameters of Vj ow _k(A = 2.1 fm -1 ) as a function of density p = 2kp/3n 2 from p /4 to 
Po- We see that all of the L = parameters, together with f ly are enhanced at lower densities. 



C. Hartree-Fock single-particle energies 

We now discuss the leading-order (Hartree-Fock) contribution to the nucleon single-particle 
energy. The second-order contributions to the quasiparticle interaction get modified through 
the resulting change in the energy-momentum relation for intermediate-state nucleons. For a 
nucleon with momentum k, the first-order (in-medium) self-energy correction reads 

k 2 v-^ - — 

tk(k F ) = — — + } (kk'sis 2 tit2\V\kk'siS2tit 2 ) 

s 2 ,t 2 ,\k'\<k F 

= ^T + ^ E( 2T + X )( 2J + X ) / ( + dpp 2 mm{2, (k 2 F - (k - 2p) 2 )/Apk} 

AM N All isjt Jmax{0,(fc-fc F )/2} 

x ( P lSJT\V\plSJT), (69) 
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FIG. 4. Single-particle energies (symbols) computed from eq. (|69p and fit (lines) with the form eq. 
(|70p characterized by an effective mass plus energy shift. 



where p = \k — k' |/2. In Fig.HJwe plot the single-particle energy as a function of the momentum 
k. In the left figure we show the results for all three NN interactions considered in the previous 
section at a Fermi momentum of kp = 1.33 fm -1 . In the right figure we consider only the 
low-momentum NN interaction with A = 2.1 fm -1 for three different densities. In all cases one 
can fit the dispersion relation with a parabolic form 

k 2 

e k = + A . (70) 

2M* K ' 

with M* the effective mass and A the depth of the single-particle potential. From the figure 
one sees that this form holds well across the relevant range of momenta k. In Fig. [5] we show the 
extracted effective mass and potential depth for the three different interactions as a function 
of the density. The energy shift A shows more sensitivity to the decimation scale A than the 
effective mass M*. At saturation density po — 0.16 fm -3 , the variation in A is about 30% while 
the spread in M* /M^ is less than 10%. Overall, the effective mass extracted from a global fit 
to the momentum dependence of the single-particle energy is in good agreement with the local 
effective mass at the Fermi surface k = kp, encoded in the Landau parameter f\. The largest 
difference occurs for the bare Idaho N 3 LO potential owing to the larger momentum range over 
which eq. ( 1701) is fit to the spectrum. 

We employ the quadratic parametrization of the single-particle energy in the second-order 
contributions to the quasiparticle interaction eqs. (!2T!) - (l24l) . This greatly simplifies the inclusion 
of the (first-order) in-medium nucleon self energy. The second-order quasiparticle interaction 
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FIG. 5. Parameterization of the single-particle energy = k 2 /2M* + A as a function of density for 
three different NN potentials. 



is effectively multiplied by the same factor M*/Mn, since the constant shift A cancels in the 
energy differences. We then compute the dimensionless Fermi liquid parameters by factoring 
out the density of states at the Fermi surface Nq = 2M*/cf/7t 2 . In Table [IV] we show the 
results at kp = 1.33 fm -1 for the Idaho N 3 LO potential as well as V£ w _ k for A = 2.1 fm^ 1 and 
A = 2.3 fm -1 . In addition, we have tabulated the theoretical values of the different nuclear 
observables that can be obtained from the Fermi liquid parameters. The quasiparticle effective 
mass of the bare N 3 LO chiral NN interaction is M*/Mm = 0.944, but this ratio increases 
beyond 1 for the low-momentum NN interactions. The inclusion of self-consistent single-particle 
energies in the second-order diagrams reduces the very large enhancement in the effective 
mass seen previously in Table IIHI Compared to the other three L = Landau parameters, 
the spin-spin interaction in nuclear matter is relatively small (Go — 0.35 — 0.58). Despite 
the strong repulsion in Fq that arises from the second-order ph diagram, we see that nuclear 
matter remains unstable against isoscalar density fluctuations (Fq < — 1), and this behavior is 
enhanced in evolving the potential to lower resolution scales. The nuclear symmetry energy 
/3 is weakly scale dependent and we find that the predicted value is within the experimental 
errors /3 = (33 ± 3) MeV. Partly due to the rather large effective mass M* /M^ at the Fermi 
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VN3LO 
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1.39 


1.59 
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0.58 


0.47 


0.96 


-148 


30.5 


0.12 


y2.3 
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-1.77 


0.54 


1.98 


2.07 


0.36 


0.74 


0.80 


0.72 


1.12 


-152 


32.5 


0.07 


t/2.1 

low— k 


-1.98 


0.58 


1.94 


2.14 


0.38 


0.83 


0.87 


0.80 


1.13 


-191 


31.8 


0.07 



TABLE IV. Sum of the first- and second-order contributions to the dimensionless Fermi liquid parame- 
ters for the Idaho N 3 LO potential and two low-momentum NN interactions VJ^_ k at kp = 1.33 fm -1 . 
Hartree-Fock self-energy insertions, as parameterized in eq. ()70p . are included in the second-order 
diagrams. 



surface, the anomalous orbital g-factor comes out too small compared to the empirical value 
of Sgi = 0.23 ± 0.03. The L = spin-isospin Landau parameter G' is quite large for the 
low-momentum NN interactions. Using the conversion factor (g n N /2Mtv) 2 = 1.9 fm 2 one gets 
the values g' NN = 0.67, 0.75, and 0.77 for the bare N 3 LO chiral NN interaction and evolved 
interactions V^_ k and V^_ k respectively. These numbers are in good agreement with values 
of 

9nn ~ 0-6 obtained by fitting properties of giant Gamow- Teller resonances. 
The above results highlight the necessity for including three-nucleon force contributions to 
the quasiparticle interaction both for the bare and evolved potentials. In fact it has been 
shown that supplementing the (low-momentum) potentials considered in this work with the 
eading chiral three-nucleon force produces a realistic equation of state for cold nuclear matter 
4(. The large additional repulsion arising in the three-nucleon Hartree-Fock contribution to 
the energy per particle should remedy the largest deficiency observed in present calculation, 
namely the large negative value of the compression modulus /C. A detailed study of the effects 
of chiral three-forces (or equivalently the density-dependent NN interactions derived therefrom 



34N36| ) on the Fermi liquid parameters is presently underway. 



IV. SUMMARY AND CONCLUSIONS 

We have performed a complete calculation up to second-order for the quasiparticle interac- 
tion in nuclear matter employing both the Idaho N 3 LO chiral NN potential as well as evolved 
low-momentum NN interactions. The numerical accuracy of our results is on the order of 
1% or better. This precision is tested using analytically-solvable (at second order) schematic 
nucleon-nucleon potentials emerging from scalar-isoscalar boson exchange and modified "pion" 
exchange. We have found that the first-order approximation to the full quasiparticle interac- 
tion exhibits a strong scale dependence in / , /g, and /{, which decreases the nuclear matter 
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incompressibility /C and increases the symmetry energy /3 and anomalous orbital g-factor 5gi 
as the resolution scale A is lowered. Our second-order calculation reveals the importance of the 
hole-hole contribution in certain channels as well as the strong effects from the particle-hole 
contribution for the fo and f\ Landau parameters. The total second-order contribution has a 
dramatic effect on the quasiparticle effective mass M*, the nuclear matter incompressibility K, 
and symmetry energy /3, as well as the Landau-Migdal parameter g' Q that governs the nuclear 
spin-isospin response. In contrast, the components of the spin-spin quasiparticle interaction 
{do,9i) are dominated by the first-order contribution. We have included also the Hartree-Fock 
contribution to the nucleon single-particle energy, which reduces the second-order diagrams by 
about 30% (as a result of the replacement — > M*). The final set of L — 0, 1 Landau 
parameters representing the quasiparticle interaction in nuclear matter provides a reasonably 
good description of the nuclear symmetry energy and spin-isospin collective modes. Our 
calculations demonstrate, however, that the second-order quasiparticle interaction, generated 
from realistic two-body forces only, still leaves the nuclear many-body system instable with 
respect to scalar-isoscalar density fluctuations. Neither the incompressibility of nuclear mat- 
ter K nor the anomalous orbital g-factor Sgi could be reproduced satisfactorily (without the 
inclusion of three- nucleon forces). A detailed study of the expected improvements in the quasi- 
particle interaction resulting from the leading-order chiral three-nucleon force is the subject of 



an upcoming investigation 



19]. 
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